What future does the universe have? 



(N 

o 
o 

(N 
-i— > 

o 
O 

m 

(N 



B. Hoeneisen 

27 March 2001 
Universidad San Francisco de Quito 
PO Box 17-12-841, Quito, Ecuador 

brucel@fnal.gov 



> 
(N 

in 
O 

(N 
O 

Of 

6 

CO 



13 



Abstract 

We discuss the future evolution of the universe in the light of recent 
observations. The apparent luminosity vs redshift of supernovae favor 
an accelerating universe. However an Einstein-de Sitter critical uni- 
verse should not be ruled out yet. 
universe: evolution 

Let us try to understand the future of the universe in the light of recent 
observations. Consider this beautiful equation: 



da 
H n dt 



h fife + 

a 



It describes the evolution of a homogeneous and isotropic universe accord- 
ing to General Relativity if the density of non-relativistic matter dominates 
radiation (and other sources of pressure). The first three terms can be ob- 
tained from Newtonian physics and are proportional to the kinetic energy, 
potential energy and total energy respectively. The last term is due to the 
cosmological constant A of General Relativity. The expansion parameter a(t) 
is normalized to a (to) — 1 today. The Hubble constant is 



H = 



1 da 
a dt 



100/io • km • s _1 Mpc _1 



he 



j , 9.778Gyr 
with ho = (0.71 ± 0.07) x Note that (0) at the present time is 

1 = Q m + fifc + Q\ 



(2) 



(3) 



so the evolution of the universe on large scales is determined by just two 
independent parameters, e.g. Vt m and VL^. is the present density of 
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non-relativistic matter (barionic and dark) in units of the critical density 
p c = 3H$/8nG. 

kr 2 
_n a 

where fc = +1, -1 or 0. If k — +1, space has the geometry of the 3- 
dimensional surface of a sphere of radius a$ x a(i) in 4 dimensions. If k — — 1, 
space has negative curvature. If = then space is flat. Finally 

Ac 2 

We consider the four universes listed in Table ^ (£l m , ^a) = (1, 0) is the 
Einstein-de Sitter critical universe with flat space that expands for ever with 
a velocity approaching zero; (1.5, 0) is overdense and has a space of finite 
volume that expands to a maximum and then collapses into a Big Crunch; 
(0.3,0) is underdense so matter can not halt the expansion; and (0.3,0.7) 
has flat space and in the future will expand exponentially for ever due to the 
cosmological constant. 

The purpose of this note is to estimate the likelihood that each of these 
universes corresponds to ours, and to understand in a simple way the issues 
involved. To obtain this estimate we compare observations with predictions 
using the "chi-square" function: 

2 \ (observation^ — theoryj 2 
z — ' error, 

i 1 

This function is written for the case when the observations are independent 
so that a judicious choice of observations is needed. In the end we shall see 
that the main source of uncertainty of this analysis is the estimate of the 
systematic errors. The contributions to the \ 2 are listed in Table |IJ 

Let us begin with observations on galaxies. We have developed a simple 
model of the hierarchical formation of galaxies 0]- The simulations were done 
with ho = 0.6. After adjusting three parameters (the amplitude A and slope 
n of the power spectrum of primordial density fluctuations, and Q m ) this 
model is in quantitative agreement with the following observations: the Tully- 
Fisher, Faber- Jackson and Samurai relations, the Schechter distribution (two 
parameters), the galaxy-galaxy correlation (two parameters), the fluctuation 
in galaxy counts, the fluctuation in the cosmic background radiation, and the 
peculiar velocities of galaxies. To compare the model with observations we 
have obtained a \ 2 with 7 terms. For each Q m we have minimized the x 2 by 
varying 2 parameters: A and n0]. Therefore for each Q m we are left with 5 
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^fc 


age 


contributions to \ 2 
G B age Q m S 


Total 

x 2 


CL 


1.0 0.0 0.0 
1.5 0.0 -0.5 
0.3 0.0 0.7 
0.3 0.7 0.0 


9.2Gyr 
8.5Gyr 
11.2Gyr 
13.3Gyr 


3.2 0.0 2.9 5.4 6.1 
3.7 13 3.8 16 12.0 
15 25 1.0 0.0 1.2 
12 0.0 0.1 0.0 0.0 


18 
49 
42 
12 


6% 
0.0% 
0.0% 

28% 


degrees of freedom: 
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Table 1: Comparison of 4 universes with observations. Density perturbations 
are assumed to be scale-invariant (n = 1): keeping n as a free parameter does 
not change the conclusions significantly, see the text. The contributions to 
the x 2 are from (G) properties of galaxies (see text); (B) Boomerang and 
Maxima observations; (age) age of the oldest globular clusters; (fl m ) density 
measurements; and (S) supernovae. (CL) is the probability that observations 
are consistent with the model. 

degrees of freedom. The corresponding x 2 's for the four universes are 3.2, 3.7, 
7.6 and 9.2, and the spectral indices n are 1.0, 0.9, 1.2 and 1.2, respectively. 
We note that observations on galaxies favor a universe with Q m « 1, ~ 
(the Einstein-de Sitter universe) and n « 1 (scale invariant spectrum) 
which is quite remarkable! In fact, for Q\ = 0, we obtain agreement with 
observations provided Q m > 0.25 and 1.1 — 0.3Q m < n < 1.4 — 0.2fi m with 
95% confidence (assuming the model is correct and the simulation volume 
(92Mpc) 3 is sufficiently large) @. So all four universes considered in Table p] 
are compatible with these observations on galaxies. However, if we restrict 
the slope n to 1 as suggested by the Boomerang and Maxima observations 
discussed below, then the x 2 ' s f° r 6 degrees of freedom are as shown in 
Table p] column G. Note that the two low density universes are disfavored by 
observations on galaxies if we set n — 1: the velocities of circular orbits of 
L* galaxies are too low, the fluctuations in galaxy counts are too low (unless 
n is tilted to 1.25), and the peculiar velocities of galaxies are too low (unless 
n is tilted to 1.3) @. 

Let us turn to the Boomerang and Maxima balloon-borne experiments. 
By observing the first acoustic peak of the fluctuations in the cosmic back- 
ground radiation, the Boomerang collaboration obtains Q m + ^a with 95% 
confidence in the range from (0.88 — 1.12) to (0.97 — 1.35) depending on the 
assumed priors and parametrizations||. We translate this into fl m + Q\ = 
1.08 ± 0.16 (all errors in this article are one standard deviation or 68% con- 
fidence level unless otherwise stated). The Boomerang collaboration also 
mentions that "with reasonable priors we find" Q m + = 1.07 ± 0.06. 
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These errors are statistical. We would like to add systematic errors, but 
these are hard to come by in cosmology A full analysis is given in ||. The 
Maxima experiment obtains Q m + Q\ = 1.0 a g 5 30 at 95% confidence ||. To be 
somewhat conservative we finally take Q m + Q\ = 1.0 ± 0.1 ± 0.1 where the 
first error is statistical and the last term is a rather arbitrary estimate of the 
systematic errors. The corresponding contributions to x 2 are listed in Table 
|l] column B. Note that the remarkable Boomerang and Maxima observations 
already rule out the two non-flat universes listed in Table [I]. Let us mention 
that the Boomerang experiment obtains a slope of the power spectrum of 
primordial density perturbations n = 1.0 ± 0.1 [Q which is consistent with 
scale invar iance. Maxima obtains n = 1.08 ± 0.10. 

Let us consider the age of the universe as inferred from the evolution of 
globular clusters. The estimates]]]] range from 18Gyr to HGyr with addi- 
tional ±10% errors of various origins. Thus current estimates are 14 ± 2Gyr, 
"with a possible systematic error of similar size" |1| . A recent study obtains 
11.5±1.3Gyr[g]. We take 14±2±2 where the first error is statistical and the 
second is a rather arbitrary estimate of the systematic error. In Table [1] we 
show the calculated ages (for the central value of H$) and their contributions 
to the x 2 - 

Now we consider measurements of the density of the universe. A compi- 
lation of methods and results can be found in The "luminosity xM/L" 
methods yield Q m in the range 0.1 to 0.4. These measurement do not in- 
clude extended dark matter halos between galaxy clusters. In our model of 
galaxy formation [0, |7j the galaxies and clusters of galaxies have halos with 
density run oc r~ 2 that extend, a grosso modo, out to the halo of the "next" 
galaxy or cluster of galaxies (once the peculiar motion is removed). This is 
possible at all times, in spite of the expansion of the universe, because of the 
ongoing hierarchical formation and merging of galaxies. For example, if clus- 
ters of galaxies have radius R, and the halos of dark matter with density run 
p oc r~ 2 extend out to r « 3-R, then Q m measured by the "luminosity xM/L" 
method is low by a factor m 3. Such extended halos are expected]^, [7|]. Is 
there any evidence to rule them out? For these reasons I believe that the 
quoted measurements of Q m are really lower bounds. The "baryon fraction 
in galaxy cluster" method obtains Vt m in the range from 0.15 to 0.35 but also 
does not include the halos of dark matter between clusters. The "large-scale 
velocity flow" methods yield results in the range 0.1 to 1[|TJ. In particular 
the study by Idit Zehavi and A. Dekel|| obtains 0.4 < Q m < 1.1 for Q\ = 
or 0.2 < Q m < 0.9 for fl\ = 0.7 at 99% confidence. The lower bound as- 
sumes ho = 0.75 and n = 1.1, while the upper bound assumes h = 0.55 and 
n = 0.9. From the evolution of the number density of x-ray clusters it is 
inferred that fl m ~ 0.74 while Q m < 0.3 is rejected with 95% confidence||. 
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For these reasons, and because the \ 2 corresponding to galaxies already in- 
cludes peculiar velocities, galaxy clustering, and velocities of circular orbits 
which are used to measure fl m , we take a conservative VL m = 0.3igJ. The 
corresponding contributions to the x 2 are presented in Table |]. 

Last, but not least, we consider the Hubble diagram of relative luminosity 
(magnitude) vs red-shift of supernovae of type la. This is a very interesting 
measurement. The results obtained by the "Berkeley Supernova Cosmology 
Project" |TI], [IT, [12], [13| have been confirmed by the "High-z Supernova Search 



Team" |14|, |TJ, [L6| . At high redshift the supernovae are fainter than expected 
(from the extrapolation of low redshift supernovae) yielding 0.75Q\ — Q m = 



0.25±0.16(stat.)±0.09(identified syst.) ±0.47(evolution)g(| |TJ in the range 
of interest. The corresponding contributions to the \ 2 are listed in Table |I] 
column S. For a flat universe this measurement translates to Q m = 1 — Qa = 
0.28i^(stat.)t[j;^(identified syst.) ± 0.27(evolution)[|Tg, [I3J. The system- 



atic error is dominated by evolution. The errors labeled "evolution" above 
correspond to an error in magnitude of ±0.2 for high redshift supernovae. 
The preliminary estimate of this error in [13] is ±0.2 mag. The correspond- 
ing estimate in |16[ is ±0.17 mag. In comparison note that the correction 
for extinction is of order 1 mag with and error « ±0.15 mag. The challenge 
is to reduce the error due to evolution, and the jury is still out. Are the 
high-redshift supernovae 10% fainter than expected because of "gray" dust, 
or evolution of supernovae, or evolution of the host galaxy, or gravitational 
lensing, or selection effects, or increased extinction in the host galaxy at high 
redshift, or is it really due to a cosmological constant? Or more mundane 
challenges that seem to be under control (but who knows?): how does the 
luminosity calibration or zero change with redshift?; how does noise sub- 
traction or the subtraction of the luminosity of the host galaxy depend on 
luminosity or redshift? Radio and optical gravitational lensing observations 
set the limit —1.78 < — < 0.27 at 95% confidence ||1 7|| . This limit cor- 



responds to —0.39 < f2y\ < 0.64 assuming a flat universe |T7|] which is already 
cutting well into the allowed region of the supernovae observations. Note 
that the size of the estimated systematic error due to supernovae evolution 
can completely change the conclusions of this study! If the error is half the 
estimate given above, then the Einstein-de Sitter universe is ruled out. 

The probabilities that observations are consistent with the model uni- 
verses is shown in column CL of Table [I]. These probabilities are propor- 
tional to the likelihood that the models are correct given the observations 
(if in advance we have no good reason to prefer one model over the other). 
We note that the two non-flat universes are excluded by several observa- 
tions. The preferred universe has (Q m ,Q^) = (0.3,0.7). Note that only one 
of the observations (almost) rules out the Einstein-de Sitter universe, and in 
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cosmology it is best to keep an open mind. 

To understand the accelerating universe (fl m , ^a) = (0.3, 0.7) let us com- 
pare it with the Einstein-de Sitter universe (1,0). This latter universe has a 
present age to — 2/3-ffo = 9.2Gyr; the expansion factor is a = [3H t/2] 2/3 ; the 
present distance to the horizon is <if/(t ) = 2c/ H = 3ct ; and the co-moving 
distance to the horizon d H {t)/a(t) will continue increasing as oc t 1 / 3 so that 
new galaxies will always be entering the horizon. The universe with a cos- 
mological constant, (0.3, 0.7), has a present age to = 0.97/H = 13.3Gyr; the 
expansion factor will grow exponentially in the future as a oc exp [v^a-^o^] > 
the present distance to the horizon is dnito) = 3Ac/H = 3.3ct ; and the 
co-moving distance to the horizon dn{t)/a(t) will approach 4.5c/ H so that 
after a few I/Hq no new galaxies will enter the horizon. In addition, the 
galaxies that are within the horizon will become exponentially redshifted 
and dimmer: their relative luminosity will decrease as oc exp [— Sy/Q^Hot] 
so that after a few 1/Hq only the galaxies of the local group, which are 
gravitationally bound, will remain visible. 

Finally let us calculate the time it will take to exhaust the fuel of the 
universe, i.e. mainly hydrogen (and also the elements from helium to man- 
ganese). We assume that hydrogen is burned at the rate we see today. The 
luminosity density of the universe is Lb = (2.0 ± 0.4) x 10 8 hoL sun Mpc~ 3 ^ 
with L sun pa 3.85 x 10 26 W. The energy released per proton is ~ 8.8MeV. The 
density of hydrogen in the universe is pu = 0.76p# pa 0.76 x 0.019p c0 //io- 
From this data we obtain a time to exhaustion of hydrogen at the present 
rate of consumption of ~ 4000Gyr, or about 300 times the present age of the 
universe. 

In summary, we have considered four universes. From Table |l] we note 
that one of them is preferred by observations: the low density universe with 
a cosmological constant. The conclusion is tentative pending a full study of 
possible evolutionary effects of supernovae and a positive detection of the 
cosmological constant by an independent method. The Einstein-de Sitter 
universe can not be ruled out yet. If the cosmological constant is not zero 
the question is: why are terms with fl m and Q\ in Equation [I] of compa- 
rable magnitude today, when one dominated by many orders of magnitude 
in the past, and the other will dominate in the future? This is the "Why 
now?" problem. Another question: why is the cosmological constant so small 
compared to predictions, yet (apparently) not zero? 

Let us assume that the universe is indeed (fl m ,Q\) = (0.3,0.7). In this 
case we live in a universe that is spatially flat and so has infinite volume; dis- 
tances between far away galaxies will expand exponentially for ever; galaxies 
will cease to enter the horizon; the ones within the horizon will dim exponen- 
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tially fast until only the local group remains visible; the galaxies of the local 
group, which are gravitationally bound, will merge into a single one; small 
perturbations 5 = (p — (p)) / (p) of the density of the universe will cease to 
grow and the hierarchical formation of galaxies will come to a halt; hydrogen 
fuel will become exhausted, stars will die leaving cinders (white dwarfs which 
will cool and stop shining, neutron stars and black holes), and darkness and 
cold will prevail for ever. No conscience will be there to know. 
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